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To the memory of A. A. Karatsuba, on 
the occasion of his 75th anniversary. 

Abstract 

If n is a positive integer, let h(n) denote the maximal value of the 
product gi<72 • • ■ Qj for all families of primes q\ < qi < . . . < qj such 
that qi + q% + . . . + qj < n. We shall give some properties of this 
function h and describe an algorithm able to compute h(n) for any n 
up to 10 35 . 

1 Introduction 
1.1 Function h(n) 

If n > 2 is an integer, let us define h(n) as the greatest product of a family 
of primes q\ < q2 < ■ ■ ■ < qj the sum of which does not exceed n. 

Let £ be the additive function such that £{p a ) = p a for p prime and 
a > 1. In other words, if the standard factorization of M into primes is 
M = q^q% 2 . . . qj j , we have l(M) =q^+q^ + ...+ q^ and £(1) = 0. If fi 
denotes the Mobius function, h(n) can also be defined as 

(1.1) h(n) = max M. 

£(M)<n 

Note that 

(1.2) £(h(n)) < n. 

From the unicity of the factorization of h(n) into primes, the maximum in 
(jl.ip is attained in only one point. It is convenient to set 

h(0) = h(l) = 1. 
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(h(n)) n >i is sequence A159685 of the OEIS (Online Encyclopedy of Integer 
Sequences). A table of the 50 first values of h{n) is given at the end of the 
paper. A larger table may be found on the authors's web sites [2j IIP], 

In [9], Landau has introduced the function g{n) as the maximal order of 
an element in the symmetric group <3 n ; he has shown that 

(1.3) q(n) = max M. 

^ ' ^ 1 £(M)<n 

The introductions of [5] and [3] recall the main properties of Landau's func- 
tion g(n) which is mentionned as entry A002809 in [12J. From (|l.ip and 
(fL3|) . it follows that 

(1.4) h(n) < g{n), (n > 0). 

In this article, we shall give some properties of h(n) and describe an 
algorithm able to calculate h(n) for any n up to 10 35 . 

1.2 Notation 

• We denote by N the set of non-negative integers. 

• The symbol p will always denote a prime number. 

• For every arithmetic function / : N — > C, we define 

(i-5) */(*)= /(p) 

p<x, p prime 

• In particular, for f(n) = 1, we will note, as usual tt(x) = iT\(x) the 
number of primes up to x. 

• For f(n) = n we define 

(1.6) ir id (x) = ^ p 

p<x, p prime 

• We denote by pj the j'-th prime and we set <tq = 0, Nq = 1 and, for 
3 > 1. 

(1.7) Oj = TTid{Pj) = Pi + P2 + • • • + Pj, Nj = pip 2 ..-Pj. 

In §[3j for all j > 1, we shall prove that h(cTj) = Nj. 

• If m is an integer, we denote by m* the smallest prime p satisfying 
p > m and, if m > 2, by *m the largest prime p satisfying p <m. 

• P + (m) (resp. P~(m)) will denote the largest (resp. smallest) prime 
factor of m > 2. It is convenient to set -P + (l) = — oo, P~(l) = +oo. 
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• oj(n) is the number of distinct prime factors of n and Q(n) the num- 
ber of prime factors of n, counted with multiplicity. u(n) is Mobius's 
function. 

• For x > 1, log 2 (x) = log log x. 

• Li is the integral logarithm defined for x > 1 by 

/•1-e rx ^ (loRx) n 
Li(x) = lim / + / — ■— = 7 + log 2 x + > v & / 
£ ^o, e >o J Q J 1+£ logt B2 ^ nn\ 

where 7 is Euler's constant. 
1.3 Functions hj(n) 

For n > 0, let k = k(n) be the non-negative integer defined by 
(1-8) a k = TT id (pk) <n< TTidipk+i) = <?k+i- 

It is the maximal number of prime factors of h(n). For < j < h = k(n), 
let us set 

(1.9) hJn) = max M 

J e(M)<n 

At(M)^0, u{M)=j 

where co(M) is the number of prime factors of M. For n > 0, we have 

(1.10) ho(n) = 1 
while, for n > 2, we have 

(1.11) h x {n)=*n>2. 
Note that 

(1.12) £(hj(n))<n. 

In §[6j we prove that, for all n's, the sequence hj(n) is increasing on j, so 
that 

(1.13) h(n) = h k (n), (n > 0). 

Our proof is not that simple. A possible reason is that this increasingness 
relies on the properties of the whole set of primes V . Let V' be a subset of 
V and M-pi the set of integers whose prime factors belong to V' . We may 
consider 

(1.14) hJn,P')= max M . 

M&N v ,,l(M)<n 
H(M)^0, uj(M)=j 
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By choosing P = {2, 3, 11, 13, 17, 19, 23, . . .} = V \ {5, 7}, we observe that 



/i 2 (24, P) = 11 • 13 = 143 > /i 3 (24, P) = 2 • 3 • 19 = 114. 

In §[H we give an upper bound for hj(n) which will be useful in §[U] 
where our proof of the increasingness of hj is given. In f j 1 . 9 [) . hj(n) can be 
considered as the solution of a problem of optimization with prime variables. 
The upper bound of hj{n) is obtained by relaxing some constraints so that 
certain variables are no longer primes, but only integers. 

1.4 Elementary computation of h(n) and hj(n) 

The naive algorithm described in |5j to compute g(n) can be easily adapted 
to calculate h{n) for 1 < n < N. Note that, for the prime factors of h(n), 
Corollary 13.11 below furnishes the upper bound 

P + (h(n)) < p k{n)+1 +P k {n)+2- 

It also can be adapted to compute hj(n). For r > j > 1 and n > o~j, let us 
define 

h^\n) = max M . 

J P+(M)<p r , £(M)<n 

M(Af)^0, u(M)=j 

We have the induction relation 

h { ; +1 \n) = nuK(fcW(n),PH-iftJ P i(n-ft4.i)). 

Indeed, either p T +i does not divide h^" +1 \n), and h^ +1 \n) = hp(n) holds, 
or Pr+i divides h^ +l \n), and h^ +l \n) = Pr+ih^^n— p r +i), which implies 
n > Pr+l + (Tj-l- 

Moreover, if p r > n, we have hj (ri) = hj(n), hr (ri) = N r and h± (n) = 

*n for n < p r while, for n > p ri h^\n) = p r holds. So, we may write 
algorithm [TJ which has been used to calculate the table in appendix. The 
merging and pruning method described in [21 §2.2] can be used to improve 
the running time. 

In § [HJ a more sophisticated algorithm to calculate h{n) is given. It is 
based on a fast method to compute Wid(x), which is explained in § 



2 Some lemmas 

Lemma 2.1. If m > 2 is an integer, let us denote by m* (resp. *m) the 
smallest (resp. largest) prime p satisfying p > m (resp. p < m). Then 

11 7 

m* < — m and *m > — m 
~ 8 ~ 10 

hold. 
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Algorithm 1 Computation of hj{n) for 2 < n < nmax and 1 < j < k(n) 
Procedure ComputeHj (nmax) 

r = 1; p = p r ; kmax = k(nmax); pmax = pkmax+i + Pkmax+2 
while p < pmax do 

for n from a r to nmax do 

H[r,n] = N r 
jmax = min(r — 1, kmax) 
for j from jmax by —1 to 2 do 

for n from nmax by —1 to p + <Jj—\ do 

H[j, n] = m&x(H[j, n],p * H[j — 1, n — p\) 
for n from p to nmax do 

F[l,n] =p; 
r = r + 1; p = p r 



Proof. We use the result of |6|: for x > 396738, the interval \x,x + nr . x ^ 1 
contains a prime number. As 396833 is prime, we deduce that, for pi > 
396833, 

, s Pi+l 1 1 11 10 

2.1 ?1±±- < i + < i + < 1.00025 < — < — . 

Pi - 25 log pi ~ 25 log 2 396833 8 7 

If m is prime, m* = *m = m holds, while, if m is not prime, we define pi 
by pi < m < pi+i; we have m* = Pi+i, *m = pi, 

m* < pj+i ^ *m > pi ^ pj 



m pi + 1 ^ m pi+i - 1 pi+i 

and, if pj > 396833, the result follows from (|2.ip . Finally, it remains to check 
that 

ll±^< l l and — >L 

Pi + l ~ 8 pi+i - 1 " 10 

hold for all pj's satisfying 2 < pi < 396833. □ 

Lemma 2.2. Let p < p' be two primes. There exists a third prime p" 
satisfying 

(2.2) p + p < p" < pp — p + 1. 

Proof. Let us show that p" = (p + p')* satisfies (|2,2p . By Lemma 12.11 it 
suffices to prove that ^(p + p') < pp' — p + 1, i.e: 

(2.3) ^( 8 _ll_I£ + A)>o. 

\ p p pp / 

If p > 3 and p' > 5, we have ^- + ^<^ + ^<8and $£5$ holds. Similarly, 
if p = 2 and p' > 11, the inequality y + ^<T + TT <8 im P lies <E3]). In 
the three remaining cases, p = 2 and p' G {3,5,7}, it is easy to check that 
p" = {p+p')* satisfies ([22]). □ 
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Lemma 2.3. Let p and p' be two prime numbers satisfying 3 < p < p' and 
pp' 7^ 15. There exists a prime p" such that 

5 

(2.4) p + p' <p" < -pp' - p. 

6 

Proof. The proof is similar to the one of the preceding lemma. From Lemma 
I2.1| to show that p" = (p + p')* satisfies (j2.4f) . it suffices to show that 
^-(p + p') < \pp' —p, i-e. 33/p + 57/p' < 20, which evidently holds for p > 3 
and p' > 7. □ 

Lemma 2.4. For all i > 2, the following inequality 

(2.5) pi +pi-i < P2i-\ 

holds. Moreover, let b be a positive integer; there exists a positive integer 
io = io(b) such that we have 

(2.6) pi + pi-i < p 2 i-b for i > i (b). 
The table below gives some values of io(b) 



b = 


1 2 


3 4 


5 


6 


7 8 


9 


10 


12 


13 


18 


30 


3675 


io = 


3 4 


7 8 


18 


19 


27 28 


36 


39 


50 


53 


85 


149 


33127 



Proof. We start from the two inequalities 

(2.7) pi < i(logi + loglogi-a), (a = 0.9484, i > 39017), 

(2.8) pi > i(log i + log log i — 1), (* > 2) 
which can be found in [7]. From (|2.7p . it follows that 

(2.9) Pi-x+Pi < (2»-l)(log* + logIogi-a), (t > 39018) 
while, if i > max(2, 6), which implies 2i — b > 2 and i>b, (|2.8p gives 

(2.10) p 2 i-b> (2i-6)(logi + log2 + log( ) +loglogi-l) 



2? 



By using the inequality log t < t — 1 , we get 



,'2i-6\ ( 2i \ ( 2i \ b 

log ^7— = "log 77-^ > - 7T-^ " 1 



2i j °\2%-bj~ \2i-b J 2i - b 

and (jZHHD yields 

(2.11) p 2i „ b > (2i-6)(Iogi + log log i + log 2- 1) - b. 



Under the condition 



(2.12) i > max(39018, b), 

the substraction of (|2.9p from (|2.11|) gives 

P2i-b - Pi-1 - Pi > (logi + log log i + log2)(l - b) 

+ 2z(log 2 - 1 + a) - log 2 - a 



(2.13) 



> (log i + log log i + log 2) 



1.283 i - 1.642 
log i + log log i + log 2 



(6-1) 



Now, the two functions t t- > t/(logt + log log t + log 2) and t i— )• — l/(logi + 
log log t + log 2) are increasing for t > e 2 ; choosing i\ = 39018 and 



(2.14) 



1 + 



1.283 h - 1.642 
log i\ + log log ii + log 2 



L3675.52 . 



3675 



shows that, for £ > i\, (|2.12p is satisfied and that in (|2.13p . the bracket is 
positive. Therefore, (|2,13p proves pi + Pi-i < P2i-3675 for i>i\= 39018. 

To determine the entries of the table, for all i's up to 39018, we have 
calculated b{ = 2i — 1 — n(pi +pj_i) which is the smallest integer such that 
Pi-1 +Pi < P2i-bi- Further, for each b in the table, we have determined io(b) 
which is the smallest integer iq such that, for io(b) < i < 39018, b{ >b holds. 

As io(l) = 3, for alH > 3, pi < pn— i holds. So, (|2.5p follows from 

= 3 + 2 = 5 =p 3 . □ 

Lemma 2.5. Under Riemann hypothesis, for all x > 41 we have 



(2.15) 



\7T id (x)-U(x 2 )\ < — X 3 / 2 log 
247T 



X. 



Proof. Let us define r{x) by ir(x) = Li(x) + r(x) and assume the Riemann 
hypothesis. Then cf. [TQ (6.18)] : 



(2.16) 



1 



r(x)| = \tt(x) — Li(x)| < — -v/xlogx (for x > 2657). 



^7T 



Let us denote xq = 2657. Then, from (|1.6p . Stieltjes's integral gives : 



Kid{x) = Tr id (x ) 



t d[7T(t)] 



7r i(i (xo)+ / td(Li(t))+ / id[r(i)] 



= n id {x )+Li{x 2 )-U{x 2 ) + tr(t)\ x Xo - I r(t)dt. 
With (f2~T6|) . it comes 



|7Tj d (x) - Li(x 2 )| < |7r id (x ) - Li(xo) - x r(x )| + 



3 / 2 log x log i 



8vr 



+ 



8tt 



dt 
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and, using J \ft\ogt = j^ 3 / 2 ^logt - ^ , 

(2.17) W id (x) -U(x 2 )\ < — x 3 / 2 logx-— X V 2 
v ' v 71 - 24vr & 18vr 



+ kid(x ) - Li(xg) - x r(x )| - -i-x 3 / 2 log^o + -^-x^ 2 . 

1Z7T 187T 

The computation of 



r(z ) = vr(x ) -Li(x ) = 384 -399.59681... = -15.59681... 
vr ld (x ) - Li(sg) = 464 653 - 480610.2863 . . . = -15957.2863 . . . 

and (|2.17|) imply for x > xq, 

5 1 5 

WiAx) - U(x 2 )\ < —x 3/2 logx x 3/2 - 740.023 .. . < — z 3/2 logx. 

' v ;| ~ 24vr & 18vr ~ 24vr B 

which proves (|2.15p for x > xq = 2657. It remains to check (|2,15p for 
41 < x < 2657; by setting 

5 5 
= Li(x 2 ) - — x 3/2 logx, f 2 (x) = U(x 2 ) + — x 3/2 logx, 

it is equivalent to check 

(2.18) h{x) < TT id (x) < f 2 (x) 

for 41 < x < 2657. One remarks that fi and f 2 are increasing for x > 2. 
Therefore, to prove (|2.18p . it suffices to check that for every prime p satisfying 
41 < p < 2657 we have fi(p') < TTid(x) < hip) where p' is the prime 
following p. □ 

Note that, in the range [2. .2657], TTid(%) — Li(x 2 ) has several changes of 
sign, the smallest one being for x = 110.35 . . . 



Lemma 2.6. Let z and u be two real numbers satisfying z > 1 and \fz < 
u < z. Suppose that we have precomputed the tables prime, piftab and pi. 
The first two tables are indexed by the integers k, < k < vr(u), and the 
third one by the integers t, < t < u. 

• prime [A;] contains p^ (po = l). 

• piftab [A;] contains irf(pk). 

• pi[i] contains vr(t). 
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Then the sum 

(2-19) E /(«)*/ (f ) 

\fz<q<u, q prime 

may be computed in O (y/z/ log z) time. 

Proof. For g > ^/i, z/q belongs to [1, y/z) . The number of primes in this 
interval is O (y/z /log z), thus the number of values of it f(z/q) is 0(y/z/ log z). 
We group the q's for which iTf(z/q) takes the same value. Algorithm [2] carries 
out this computation. 

Algorithm 2 : Computation of the sum (|235D in O (y/z/ log z) time 

5 = 0; imin = 1 + pi[Lv^J] 
while imin < pi[u] do 

q = prime [imin] 

s = pi[z/q] 

imax = min(pi(z/ prime [s]), pi [u]) 
5 = 5+ (piftab [imax] — piftab [imin — 1]) * piftabfs] 
imin = imax + 1 
return 5 



Let us give some words to convince of the correctness of algorithm 2 : let 
us note s = %(z/q). Then p s is the largest prime < z/q. For q' prime, q' > q, 
we have -Kf(z/q') = 7Tf(p s ) = piftab[s] if and only if z/q' > p s i.e. q' < z/p s , 
in other terms, n(q / ) < ir(z/p s ). Thus the largest prime q' in the range [q..u] 
such that iTf(z/q') = iTf(p s ) is pi where i = min(7r(z/p s ), ir(u)). □ 

3 First results 

Proposition 3.1. Let j be a positive integer and o~j and Na be defined by 
([177]) . We have 

h(a j ) = Nj. 

Proof. It is easy to see that h(o~\) = h(2) = 2 = N± and h(o~2) = h(5) = 6 = 
N2- Now, we may suppose that j > 3, i.e. pj > 5 and we set p = pj/logpj. 
The function t \— > tj log t is increasing for t > e and, since 2/ log 2 < 5/ log 5, 
we have, for 1 < i < j, pi/logpi < p and for i > j, pi/logpi > p; in other 
words, i — j and pi/logpi — p have the same sign. 

Let M be a product of r distinct primes, M = Q1Q2 ■ ■ - Q r , with r > 0. 
After a possible simplification by s primes (0 < s < min(j, r)), we may write 
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with u = r — s, v = j — s and 

Pkl < Pk 2 < ■ ■ ■ < Pk v < Pj < Ph <Pj 2 < ■■■ < Pj u - 

Let f(M) = i(M) — plogM. From the definition of £, the function / is 
additive and we have 

U V 

(3.1) /(M) - f(Nj) = Y,(jPk ~ PlogftJ " I>* " P l °ZP^ ^ 

i=i i=i 

since each term of the first sum is non-negative while, in the second sum, 
each term is non-positive. 

From (jl.ip . since £(Nj) = aj, in order to prove that h((Tj) = Nj, we must 
show that, for all squarefree number M satisfying £(M) < Oj = l(Nj), we 
have M < Nj. But, for such an M, (^TT) yields 

f(M) = £(M) - plogM > f(Nj) = £{Nj) - plogNj = a, -plogNj 

whence 

M f£(M)-a i \ 

which completes the proof of Proposition 13.11 □ 

Proposition 3.2. Let r and j be two positive integers and Oj, Nj and hj be 
defined by (flTT)) and ([L9]) . We have 

(3.2) hj(a j+r - a r ) = N j+r /N r = p r+1 p r+2 . ..p r +j- 
Moreover, if n > CFj +r — a r we have 

(3.3) £(hj(n)) > o-j+r — cr r . 

Proof. The proof is similar to the one of Proposition 13.11 Let us set 

(3.4) p = -^±r — an d p > = p\ogp r - p r . 

log(p j+r /p r ) 

Since, for t ^ 1, (t — l)/t < logt < t — 1 holds, we have Pj+ r > p > p r > 2. 
For a squarefree number M, we consider the additive function 

f(M) = £(M) - plogM + p'oj(M) = f(P) = I> " Plogp + P). 

p\M p\M 



We will prove that / attains its minimum in N = Nj +r /N r . From ()3.4|) . 
it follows that f(pj+ r ) = f{p r ) = and the study of the function t i— > 
i — p log t + p' shows that 

> for p < p r or p > pj +r 

/(p) < < for p r <p < p j+r 

^ = for p = p r or p = Pj+r- 



10 



Therefore, we have 



(3.5) f(M)-f{N)= Yl f(P)~ E f(P)>°- 

p\M p\M 

p < p r Or p > p j+r p r < p < Pj + r 

From (|1.9p . we have to show that, for any squarefree integer M satisfying 
£(M) < a j+r - a r = £{N) and u(M) = j = u(N), we have M < N. For 
such an M, ([3.5p gives 

£(M) - plogM + p'w(M) > £(N) - plog N + p'u(N) 

yielding 

M . f £(M)-£(N) \ < i 



iV \ p J 



which, together with £(N) = <Jj+ r — &ri proves (|3.2p . 

To prove (|3.3p . first, from (|1.9p . we observe that hj(n) > N = Nj +r /N r . 
Setting M = hj(n) in (|3,5p and noting that w(M) = u(N) = j, we see that 

M 

£{M) > £{N) + plog — > £(N) = a j+r - a r 



which proves (|3,3p , □ 

Proposition 3.3. Let n > 2 be an integer and p < p' two prime numbers 
which do not divide h(n). Then the largest prime divisor P + (h{n)) of h(n) 
satisfies 

P + (h(n)) <p + p. 

Proof. Let us assume that the set of prime factors of h(n) not smaller than 
p + p' is not empty and let q > p + p' be its smallest element. 



• If q < pp' , by setting M = *%-h(ri), we have by (|1.2p 

£{M) =p + p -q + £{h{n)) < £(h(n)) < n 
and thus, from (jl.ip . 

(3.6) h(n) > M =^-h(n), 

in contradiction with q < pp'. 

• If q > pp' , i.e. g > pp' + 1, by Lemma I2.2[ the interval [p + p', q — p] 
contains a prime; thus the prime p" = *(q — p) satisfies p + p' < p" < 
q — p < q and, from the definition of q, p" does not divide h(n). By 
Lemma |2.1[ p" = *{q — p) > jq(q — p) holds, whence 

10 pp" [10 7 

q<^rp +p=-=- — h-^ 

7 7 \ p p 
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We have p > 2, p' > 3 and p" > p+p' > 5, so that y + ^7 < t + I < 7 > 
yielding q < pp" . By considering M = ^—h(n), as in (|3.6p . we get 

h(n) >M = ^-h{n) > h(n), 

a contradiction. 

□ 

Corollary 3.1. If k = k(n) is defined by (j 1 . 8 [) - the largest prime factor of 
h(n) satisfies 

P + (h(n)) < p k+1 +p k+ 2- 

Proof. The number of prime factors of h(n) does not exceed k, so that, 
among pi , p%, ■ ■ ■ , Pk+2 there are certainly two prime numbers p and p' not 
dividing h{n). By applying Proposition 13. 3[ we get P + (h{n)) < p + p' < 

Pk+l+Pk+2- □ 

Proposition 3.4. Let n > 5 be an integer, k > 2 be defined by ([1.8p and 

j an integer satisfying 2 < j < k. Let us supose that there exists two 
prime numbers, p,p' not dividing hj-i(n), and satisfying 3 < p < p' and 
P + {hj-i{n)) > p + p' where P + (hj-i{n)) is the largest prime divisor of 
hj-i(n). Then the inequality 

6 

hj(n) > -hj-i{n) 
o 

holds. 

Proof. Let us consider two cases : 

Case 1 : pp' > 15. Let us denote by q < P + (hj-\(n)) the smallest prime 
dividing hj-i(n) and satisfying p + p' < q. 

• If q < —pp' ', we set M = — hj^\{n)\ we have uj{M) = j and £(M) = 

p + p' -q + £(hj-i(n)) < £(hj-i(n)) so that, from (fTTD . i(M) < n 
holds and (|1.9|) yields 



6 

(3.7) hj(n) >M> -hj-iin) 

o 



as required. 

if q > ^pp', 

o 

holds, and, from the definition of q, p" does not divide hj—i(n). 



5 

If q > —pp' , we set p" = *(q—p); from Lemma [2.3| p+p' < p" < q—p < q 
6 
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By Lemma |2.1| we get p" = *(q — p) > jfi(q — p), which implies 

10 pp" (10 7 

7 7 \ p p 

But p > 3, p' > 7, p" > p + p' > 10, thus p" > 11, and ^ + ^ < 

T + IT < T> yielding 9 < % setting M = '^-hj^i(n), as in 

(I3~71) . we get hj(n) >M> §/&,•_! (n). 



Case 2 : p = 3, p' = 5. 

• If P + (hj-i(n)) < 13, which implies n < 7Tjd(13) = 41, examining the 
table of Fig. 1 shows that, for n < 41, we have hj(n) > |/ij_i(n) with 
equality if and only if hj—\{n) = 5, 35, 385 or 5005. 

• If > 17, and 11 does not divide P + (/ij_i(n)), then we 

apply case 1 with p = 3, p' = 11, while, if 11 divides P + (h,_i(n)), 

3-5 6 

hj(n) > hj-i(n) > —hj-i(n) holds. 

11 5 

□ 



4 Bounding hj(n) 

Proposition 4.1. Let j > 1 and n > 0"j (where Oj has been introduced in 
(|l,7p J 6e two integers; we define r >0 by 

(4.1) Oj+ r — o> < n < (jj+ r -|-i — oy+i 
and n' fry 

(4.2) < n' = n - (a j+r - oy) < Pj +r +i - Pr+i- 
Then we have 

Pj+r+l _ Nj+ r +l 



(4.3) hj (n) < Pr+lPr+2 ■ ■ ■ Pr+j 



Pj+r+l - n' N r (p j + r+l - n') 



Proof. From its definition (|1.9p . hj(n) is a product of j primes. Let us denote 
by Ax, A2, ■ ■ ■ , A u (with < u < j) its prime factors exceeding pj+ r +i and 
by Bi, B2, ■ ■ ■ B r +i+ u the primes < pj+ r +i and not dividing hj{n); we have 

(4.4) hj (n) = N ^+iMA,...A u 

(where the product A\A% . . . A u should be replaced by 1 when u = 0) and 

(4.5) 2<Bx< ... < B r +l+ u < Pj+r+l < Pj+r+2 <Ai<...< A u . 
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Further, let us introduce v = £(hj(n)); by (|l,12p and (|4.ip . we have 

(4.6) v < n < (jj +r+ i — (J r +i 
and it follows from Proposition 13,21 f)3.3[) . that 

(4.7) v = l(hj(nj) > (jj+ r — a r . 
Moreover, (j4.4j) implies 

u r 

(4.8) v = £(hj(n) = a j+r+ i -a r + y^( A i ~ B r+i+i) ~ ^^{Bi-Pi) - B r 

i=l i=l 

Now, we consider the optimization problem (where u, r, u, A\, A 2 , ■ ■ ■ 
are fixed) 

(4.9) M = max — — 

zev f(Z) 

where T> is a subset of N r+1+it , Z_ = {Z\, Z2, ■ ■ ■ , Z r+ i +u ), 

f(Zj = Z1Z2 . . . Z r+ i +u 
and the set T> is defined by 

(4.10) Zi>pi, (l<i<r + l), 

(4.11) Zi < Z r+1 , (1 < * < r), 

(4.12) Z r+1 < Z r+l+i <Ai, (1 < % < u) 
and 

(4.13) U(Z) - R{Z) - Z r+1 + a j+r+1 -a r = v 
with 

u r 

(4.14) U(Z) = ^(A t - Z r+1+i ), R{Z) = ]T(Z, - Pi ). 

i=i i=i 

Note that, from (03]) and (j4~8]) . BeDso that ([Oil implies 

hj(n) A\A%...A U . „ AiA 2 ...A u 
4.15 = 2 <M = max 2 

^j+r+i /(#) ^e© /(£) 

IfZeV, from (OHl) . flUED and (051> . it follows that 
2=pi<Zj<,4 u , l<i<r + l + n 
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so that f(Z_) does not vanish on T> and T> is finite. Therefore, the maximum 
M. defined by (|4.9p is finite; let C be a point in T> where the maximum A4 
is attained. We shall prove that 

(4.16) U(C) = R(C) = 0. 

For that, first we claim that one of the two numbers U(C_), R(C_) vanishes. 
Indeed, assume that U(C) / and R(C) ^ 0. From $~T0\i , we have R(C) > 
0; thus there exists iq, 1 < iq < r, such that 

(4.17) C io > Pio + l>p io . 

Similarly, from (|4.12p . we have U(C_) > 0, and there exists i%, 1 < i\ < u 
such that 



(4.18) A h > A h - 1 > C r+1+h . 
Let us define C' £ N r+1+U by 

c i = c io - !> Cf'+i+n = Cr+i+n + 1, C 1 • = Ci for i ^ i , r + 1 + %\. 

To prove that (7 G £>, we have to check that C| > Pi (which follows from 
H^D ), that C' iQ < C r+1 (which follows from C[ o = C io - 1 < C r+ i - 1), 
that CV+i < C^ +1+il (since C r+1 < C r+ 1+il and C r +n-ii = C£+i-Hi - 1), 
that C' r+1+h < A h (which follows from KTfy ) and that f7(C') - R(C') = 
U{C)-R(C) (which follows from U{C) = U(C) - 1 and R(Q!) = R(C) - 1). 
Further, we have 

f(C!) C' io C' r+1+ii ^ (Cj - l)(C r+ i + j 1 + 1) 

f(C) C l0 C 7 

(4.19) = 1 _(^M <1 

since, from the definition of T> (cf. (|4.11|) and (|4.12|) ). Cj < C r +i < CV+i+ij 
holds. But (|4.19p contradicts the fact that the maximum in (|4.9p is attained 
in C. 

Let us show now that it is impossible to have simultaneously U (C_) > 
and R(C_) = 0; indeed, let us assume that U(C_) > 1 and R(C_) = (which 
implies r = or Ci = pi for 1 < i < r). We define i\ as in (|4.18p . Since 
C<ET>,we get from (^T3l) 

CV+i = (Tj+r+i - a r - V + U(C) = ((Tj+r+l - (J r +1 ~ v) + Pr+1 + U (C) 
which, by (HSJ) and U(C) > 1, yield 

(4.20) C r+ i > p r+ i + C/(C) > Pr+1 + 1. 
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We define C G N r+1+U by 

C' r+1 = C r+1 - 1, C' r+1+h = C r+ i +il + 1, C[ = Ci for ij= r + 1, r + 1 + i\. 

To prove that C' G £>, we have to check that C[ > pi for 1 < % < r + 1 
(which follows from C[ = Ci = pi if i < r and from ([4.2Up if i = r + 1), 
that C[ < C r+1 for 1 < i < r (which follows from C[ = Ci = pi < p r and 
from C' r+l = C r+ i - 1 > p r+ x, via dCTD ). that C' r+l < C' r+l+i for 1 < i < u 
(which follows from C' r+l < C r+ \ and C' r+1+i > C r+ i +i ), that C' r+1+il < A ix 
(which follows from ($JE§) and that U(C') - C' r+1 = U(C) - C r+1 ) (which 
is easy). As in (|4,19p . we have f{C_) < f(C_), contradicting the fact that the 
maximum in (|4.9p is attained in C_. 

To prove (|4,16p . it remains to show that we cannot have R(C_) > and 
U (C_) = 0. Let us suppose that R(C_) > 1 and U (C_) = 0, which implies 
u = or, for 1 < i < u, 

(4.21) C r+1+i = Ai> Pj+r+2 > Pj+r+1 + 2, 

with the help of (|43j) . From (|4.13p and (f4~71l . this time we get 

C r+ i = crj +r+ i — a r — v — R(C) < dj +r+ i - a r - (o-j +r - a r ) - R(C) 

(4.22) = p j+r+1 - R{C) < p j+r+1 - 1. 

Here we choose iq as in (I4.17P and set 

c 'i = C io ~ c 'r+i = c r+i + 1) c 'i = Ci for i^i ,r + l. 

To prove that C' G T>, we have to check that C[ > pi (which follows 
from (|4.17D ). that C' r+1 > p r +\ (which follows from C' r+l = C r +i + 1 and 
Cr+i > Pr+i), that, for 1 < i < r, C[ < C^ +1 (which follows from C[ < Ci 
and C' r+1 > C r +i), that, for 1 < i < u, C' r+1 < C r -^i-^i 1 = Ai x (which follows 
from pH and (IQ51) ) and that ^(CO+C^! = i?(C)+C r+ i (which is easy). 
As precedingly in (|4.19p . we observe that /(C') < f(C_), contradicting the 
fact that the minimum is attained in C_. 

In conclusion, we have proved (|4,16p so that Ci = pi for 1 < i < r and 
C r +i+i = Ai for 1 < i < u. Moreover, (|4.6p yields v < n, and, from (|4.13p 
and (|4.2p . we get 

Cr+1 = ^j'+r+l — °"r — V > (Tj+ r +i — Cr r — 71 = ^j'+r+l — 1 ■ 

Therefore, the maximum in (|4.9p satisfies 

^1^2 • • • A u 1 
A4 = — < 



PlP2 ■ . .p r C T+ iAiA 2 ...A u p X P2 ■ ..p r (pj+r+l ~ U 1 ) 

which, via ()4.15p . proves (|4.3p . □ 
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Proposition 4.2. With the notation of Proposition [7T7] , we have 

(4.23) h,(n) > - N > +r+l - = ^±1 
^ ' 3> ~ N r {p 3+r+l -n'Y qN r 

where q = (pj+ r +i — n')* is the smallest prime satisfying q > pj +r+ \ — n' . 

Proof. From (|4,2p . we have p r +\ < Pj+r+i — n' < pj +r+ i which implies 
Pr+i < q < Pj+r+i so that M = Nj +r+ \/(qN r ) is an integer with exactly j 
prime factors. Further, by (|4.2p . we have 

£(M) = (Tj +r+ i - a r - q < (Tj+r+x - a r - (pj +r +i - n') = n 

and, by ([O]) . %(n) > M holds. □ 

Corollary 4.1. We keep the notation of Proposition ^. 1\ if q = pj +r+ \ — n' 
is prime then 

(4.24) hj(n) = hj(a j+r+1 - a r - q) = ■ 

ql\ r 

Proof. Corollary 14.11 follows from Propositions 14.11 and 14.21 □ 

5 A parity phenomenon 

Proposition 5.1. Let k > 2 be an integer and a be an even number satisfying 
4 < a < Pk+i and hk defined by (|1.9p . We have 

(5.1) h k (cr k+ x - a) = h k (a k+1 - a- 1). 
Proof. Since n \— > h k {n) in non-decreasing, we have 

(5.2) h k (a k+ i - a) > h k (a k+1 - a- 1). 

Let us set n = o~ k +i ~ a an d note that n satisfies a k < n < a k +i so that, 
from (|1.8p . k = k(n) = k{n — 1). Let M be a positive squarefree integer 
such that i(M) < n and uj{M) = k. Such a number M is even; if not, we 
would have £(M) > 3 + 5 + . . . + p k +i = Cfc+i — 2 in contradiction with 
£(M) <n = a k+ \ — a < o~ k+ i — 4. Therefore, £(M) is the sum of 2 and k — 1 
odd numbers, so that £(M) = a k = o~ k +i + 1 = Cfc+i — a — 1 (mod 2). So, 
£(M) cannot be equal to a k+ \ — a and £{M) < a k+ \ — a — 1 holds. Thus, 
from ([1.9p . we get h k {a k+ \ — a) < h k (a k +i — a — 1), which, with (|5.2p . proves 

(EH). □ 



Proposition 5.2. Lei 6e an integer, k > 2, and q a prime number satis- 
fying 3 < q < p k . By setting m = cr k +i — q — 1, we have 

(5.3) h k -i(m) = h k _ x (m - 1) = %_ 1 (cr fe+1 - q- 2) = -|^- 
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Proof. From the table of Figure [TJ we have /ii(6) = hi(5) = 5, /12(H) = 
7i2(10) = 21 and ^(13) = 7*2(12) = 35 so that the proposition is true for 
k = 2, q = 3 and for k = 3 and q = 3 or 5. So, from now on, we assume 
k > 4. Corollary 14,11 with j = k — 1, r = 1 implies %-l( m — 1) = N k+ i/{2q) 
and, since n 1— >• h k -i(ji) is non-decreasing, it follows that 

(5.4) Ttfc-i(m) > h k -i(m - 1) = 

Let M be a positive squarefree integer satisfying ^(M) < m and w(M) = 
A; — 1. In view of (|5.4p and (|1.9p . to prove that h k -i{rn) = N k +i/(2q), it 
suffices to show that 

(5.5) M<^±1- 
V ; - 2q 

If M is odd, £(M) is the sum of k — 1 odd numbers, which implies 
£(M) = o~k = o~k+i — q = m + I (mod 2). 

So, £{M) cannot be equal to m; since, by (jl.9p . £(M) < m holds, we should 
have £{M) < m — 1; therefore, from (|1.9p . we get 

(5.6) M</i fc _i(m-l) = ^±i- 

If M is even, we have u(M/2) = A; - 2 and £(M/2) < m - 2, so that 

(5.7) M < 2hk-2{m — 2). 

• If g > 11, since we have assumed A; > 4, i.e. Pfc+i > 11, we have 

o"fc-5<m-2 = CT fc+ i - q - 3 < cr fc+ i - 14 < cr fe+ i - 10. 

By Proposition ^. ll with 7 = A;— 2, r = 2, n = m—2, n' = n — (& k — 0*2) = 
Pk+i -q + 2,we get 

h k - 2 [m -2) < 



6(q - 2) 
which, by (|5.7p . gives 

(5.8) M <J^H = ^±i^<^ii<^±i. 
V ; ~ 3{q - 2) 3q q - 2 ~ 3q 9 2q 

If q £ {3, 5, 7}, since k > 4 and Pfc+2 > P6 = 13, we have 

a k+ i - 10 < m - 2 = o-fc+i - g - 3 < cr fc+ i - 6 < 0^+2 - 17. 
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and Proposition 14.11 with j = k — 2, r = 3, n = m — 2 and n' 
n - (a k+1 - 10) = 7 - q yields 

(5.9) h k _ 2 (m-2)< A/ •-- 



30(p fc+2 +9-7) 
Using q < 7, k > 4 and Pfc + 2 > P6 = 13 gives 

_ Nfc+l p fc+2 iVfc+l 13 



30(p fc+2 +q-7) 30 p fc+2 + g - 7 ~ 30 q + 6 

13iV fc +i g < 13iV fc +i 7 iVfc+i 
30g g + 6 ~ 30<j 13 Aq 

which, together with (|5.7|) and (|5.9|) . proves 

(5.10) M<^±i- 

2g 

Inequalities (|5.6p . (|5.8|) and (|5.10p prove that (|5.5|) holds, which, with 
(|5.4p . completes the proof of f|5 . 3[> . 



□ 

Proposition 5.3. Let k be a positive integer and m = a k+ \ ~ 1/ we /lave 

(5.11) ^fc(m) = /ife(m - 1) = h k (a k +i - 2) = — |tL 

Proof. It is the same proof than for Proposition [5.21 By Proposition [3~2l with 
j = k and r = 1, we have 

fr fc (m) > h h (m - 1) = —7^- 

Further, let M be a positive integer satisfying £(M) < m and uj(M) = k. If 
M is odd, by the parity phenomenon, we have £(M) = m — 1 (mod 2) so 
that £{M) < m - 1 and M < h k (m - 1) = If M is even, we have 

M < 2/i / t_i(cjfc + i — 3) and, if fc > 2, i.e. > 7, Proposition 14.11 with 

j = k — 1, r = 2, n = tr^j — 3, n' = 2, yields 

M < 2 N k+2 = N k+1 p k+2 < 7iVfc + i iVfc+i 
" 6(p fe+2 - 2) 3 p fc+2 - 2 " 15 2 

If k = 1, it is easy to check that (|5.1ip still holds. □ 
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6 The increasingness of hj(n) on j 

Theorem 6.1. Let n > 2 be an integer and k = k{n) be defined by (jl.8p ; 
for j satisfying 1 < j < k, we have 

(6.1) fti-i(n) < -hj(n) 

and ()6.ip is an equality if and only if j = k(n) > 2 and n = o~j+i — 4 or 
n = o-j + i — 5. 

Proof If j = 1, it follows from (fTTUj) and (fl~TTj) that h Q (n) = 1, hi(n) = 
*n > 2 and hQ{n)/hi{n) < 1/2 < 5/6. which proves (|6.ip . So, from now on, 
we assume j > 2. 

The sequence (crj+r — o"r)r>o is increasing and goes to infinity. So, we 
may define rj > and n'j by 

(6.2) Vj+rj - o- Tj <n< Gj +Tj+ x - a rj +l 
and 

(6.3) n'j = n - {o-j +rj - a T] ). 

We shall consider four cases : rj < j — 4, rj > j + 3, j — 3 < Tj < j ' + 2 and 
j > 25, j - 3 < rj < j + 2 and j < 24. 

First case : rj < j — 4 

From (|2.5p and our hypothesis j > r,- + 4, we deduce 

(6-4) p r3+ i + p rj+2 < P2r 3 +3 < Pj+r, < Pj+rj+l 

and 

(6-5) p Tj+2 + Pr 3 +3 < P2^+5 < 

Let us set From (|6.2p . we get 
(6.6) < n'j = n - {a j+r] - a rj ) < Pj+^+i - p rj +i 

and applying Proposition 14.21 yield 

( 67 ) h .( n ) > ^i±J+l = N i+r ]+ l 

[7) QN rj (p j+rj+1 -n'j)*N rj - 

In view of bounding hj-i(n), we have to determine rj_i such that 

(6.8) aj-i +Tj _ 1 - a Tj _ x < n < o-j +rj _ 1 - a rj _ 1+ \. 

We shall distinguish two sub cases. 
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Sub case one, r 3 -_i = r 3 ; + 1 
Let us asume that 

(6.9) aj +r . - a rj < n < (Tj+ rj +i - cr rj+2 . 
i.e. from (|6.3p . 

(6.10) < rij- = n - (a j+r . - a r .) < Pj+ rj +i - Vr 3 +i - p rj +2- 

Note that, from (|6.4p . the right hand side of (|6.10p is positive. Then, we 
have Tj—i = rj + 1 since, from (|6.9p . 

a {j-l)+{r j +X) ~ Vrj+l = CTj+rj ~ ^r 3 +l < Vj+rj ~ ^ < « 

< (Tj+rj+1 - = 0-(j-l)+(r J+ l) + l - C7( r . j+1)+1 

holds. Via (|6,3p . this implies that 

n j _ 1 =n - {Oj-i +rj _ 1 - 0V f _i) = n - a j+rj + cr rj+ i = nj + p rj+ i. 
Applying Proposition 14.11 and noting that j — 1 + fj-i = i + f*j yield 



N rj-i (Pj-l+rj-i+1 ~ n 'j-l) N rj+1 (p j+rj+1 - n'j - p Tj+ i) 
By using (|6.7p . we get 



/ \* 



(6-H) 4^ < 



hj-i(n) (Pj+r j+ l -n'j) 



hj (n) p r . +1 (pj+r J+ i - rij - p rj+ i) 

From (|6.6p . we have pj +rj+ \ — n'- > p rj +\ > Pi = 2, so that we may apply 
Lemma [2 . 1 1 which . with the help of (|6.1ip and (|6.10p . yields 

hj-l(n) , 11 / Prj+ i \ 11 / 1 | 1 



hj(n) 8p rj +i \ p j+rj+1 - n'j - p rj+1 J 8 \p rj +i p rj +2 
If r,- > 1, ^ (■=£— + =^-) < t (s + s) < B' which P roves »• 



It remains to consider the case r, = 0, which, from (j6.9|) and (|l,8p . implies 
Uj <n < <Tj+i and /c(n) = j. 

()6.10p becomes < n', ; = n— crj < Pj+i — 5 and, by setting a = Pj+i — n'- , 
we get 

(6.12) 5 < a = - n^- = p j+ i + aj - n = <jj + i - n < pj + i 
while (|6.11|) yields 

(6.13) 3 \' < 



hj(n) ~ 2 (a - 2) 
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„ T . . 11 , , . 11 a 1112 5 

By Lemma 2.1 a < — a holds, and, for a > 12, < < — , 

y — ' ~ 8 ' ~ ' 16a- 2 ~ 16 10 6 
which, via (|6,13p . proves (|6.ip . 

Since, from (|6.12p . a > 5, it remains to study the cases 6 < a < 11. If 

a* 5 

a = 7, 9, 10, 11, it is easy to check that — r < — • 

2(a — 2) 6 

If a = 6 or a = 8, by Proposition 15.11 <j6. 12|) and (j4.23|) we have 
hj(n) = hj(<jj + i — a) = hj(aj + i — a — 1) = hj(n — 1) > 



while, by Proposition 15.21 since a — 1 is prime, we get 



hj-\{n) = hj-i(n — 1) 



N j+ i 



2(o- 1) 
yielding 

hj-i(n) _ hj-i(n - 1) (a + 1)* _ J 7/10 if a = 6 
/ij-(n) ~ /ij(n-l) ~ 2(a-l) ~ [11/14 if a = 8 

and, in both cases, — ; — — — — < — holds, which proves (16. ip . 
Aij(n) 6 

Sub case two, r^-i = r j + 2 

Now, we asume that (|6.2p holds but not (|6.9p ; thus we have 

(6.14) Oj +rj+ i - (Trj+2 <U< Oj +rj + i - a Tj +l 

and, from (|6.5|) . 

(6.15) Pi+r-,+1 - Pr^+l - Pr,+2 < < Pj+ rj +l ~ P rj +1- 
Here, we get rj-i = + 2, since we have 

cr j-l+r J +2 — 0>j+2 < < Oj+rj+l ~~ < cr j+f'j+2 — &rj+3 

by observing that, from (|6.5p . 

0"j+r j +2 - ~ (vj+rj+l ~ °>j+l) = Pj+rj+2 ~ Prj+2 ~ Prj+3 > 

holds. Now, we have 

< n' j _ l = n- (cr J _i +f . j _ 1 - cr^.J =n- <Jj+ rj +i + cr rj+2 

= n'j + <?j+r 3 - o Tj - (Tj +rj+ i + 0Yj+2 from (|6.6p 

(6.16) = n'j -pj +r . + i +Pr J+ i + < P rj +2 from (|6.15p 
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and applying Proposition 14.11 gives 



while Proposition 14.21 yields 



We set a = Pj+ rj +i - n'- and A = Pj +rj +2 ~ Pr 3 +i ~ Prj+2 so that f|6.16j) 
allows to write pj +rj+ 2 — n j-i = ^ + a > ana - we have 

( 6 17) hj-iM < Pj+r,-+2 a* 

hj(n) ~ p rj+ xp r . +2 A + a 

(|6.16p can be rewritten as 

(6.18) 2 = pi < p^+i < a = Pj+rj+i - n'j < p Tj +i + p T] +2- 

Lemma 12.11 implies a* < lla/8 and, by (|6.4p . A > holds, so that the 
homographic function 1 1— > t/(A + t) is increasing. From (|6.18p . we thus have 



< 



A + a A + + p rj+2 Pj+rj+2 

Therefore, we get 

hj-i(n) ^ UPrj+l + Pr s +2 11 ( 1 1 



/ij(n) 8 p r . +1 p r . +2 8 \p r ,+i p r , +2 



5 

which is smaller than — if r -, > 1. 

6 J ~ 



It remains to consider the case Tj = 0, r 3 -_i = 2. Formula (|6.17p becomes 
(6.19) J \ ; < ^ + 



/ij(n) 6 p J+2 -5 + a 
while ()6.18p via (|6.3p becomes 

(6.20) 2 < a = Pj+i — n'j = (Jj+i — n < 5. 

Since j ' > 2 holds, note that (|6.20p implies Uj < <7j + i — 5 < n < Oj+i — 2, 
which shows from (jl.8p that fc(n) = j. 

• If a = 5, since n = dj + i—a = Oj+\— 5, by Corollary 14. 1 1 with r = Vj = 
and g = 5, we get hj(n) = Nj+i/5, while Proposition 13.21 gives 

hj-i(a j+1 - 5) = hj-^aj+i - a 2 ) = = — - — 

Therefore, hj^i(aj + i — 5)/hj(crj+i — 5) = 5/6. 
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• If a = 4, by Proposition 15. 1| we get hj{aj+x — 4) = hj{aj+x — 5) and, 
by Proposition 15. 2| hj-xi&j-tx — 4) = hj—i(aj+i — 5) so that 

hj-i(a j+1 - 4) _ - 5) _ 5 

hj(a j+ i - 4) hj(a j+1 - 5) 6 

• If a = 3, Formula (|6.19p becomes 

hj-i(n) < Pj+2 < _7_ 5 
/»j(n) " 2(p i+2 -2) - 10 6 

since > Px = 7. 
Second case : rj > j + 3 

From flfH]) . we deduce n > (Tj +r . - <J rj+ x = o"(j-i)+(r 3 +i) - cr( rj +i) an d 
Proposition 13.21 (|3.3p . implies 

(6.21) l(hj-x(n)) > a j+rj - a r .+\. 
Let us now show that 

(6.22) q^P+(h^x(n))>p j+rr 

Indeed, if q < pj +r ._x holds, since hj-\{n) has j — 1 prime factors, we should 
have 

e(hj-x(n)) < Pj +rj -x+Pj+ rj -2 + ...+p rj +x 

< Pj+rj + Pj+rj-1 + • • • + Prj+2 = Cj+Tj — °V 3 +1 

which would contradict (|6.2ip . 

Further, among the j + 1 primes p2 = 3,p3, . . . ,Pj+2, there are certainly 
two primes p and p' not dividing /i,_i(n) and satisfying 3 < p < p' < Pj+2- 
By Lemma El p3]) . and (^22|) . we get 

(6.23) p + p < p j+1 +Pj+2 < P2j+3 < Pj+r, <q = P + {h j -i{n)) 

h-_i(n) 5 

and, applying Proposition 13.41 proves \ < — • 

hj(n) 6 

Third case : j — 3 < Tj < j + 2 and j > 25 

The proof is the same than for the second case; only, in (|6.23p . instead of 
([ZgD , we use (USD with b = 7, i = j + 2 > 27 : 

P + P < Pj+X +Pj+2 < P2j-3 < Pj+r, <q = P + (hj-x(n)). 
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Fourth case : j — 3 < r j < j + 2 and j < 24 

Here, we have rj < j + 2 < 26 and, from (|6.2p . we get 

n < <jj +rj+ i < ci5i = 5350. 

So, for < 50, (Tk < n < at+i an d 1 < j < A;, we have computed hj(n) with 
the algorithm described in Section \1. 41 and we have checked that, for j > 2, 

hj-i(n) 5 
~~ 6 

always holds, with equality if and only if j = k(n) and n = <x,+i — 4 or 
n = CTj+i — 5. □ 

Corollary 6.1. For all non-negative integer n>2, we have 

(6.24) h(n) = h k {n) 

where k = k{n) is defined by (|1.8p . 

Proof. From (jl.ip and (jl.9p we have 

h(n) = max hAn) 
0<j<k(n) J 

and Theorem 16.11 yields ho(n) < h\(n) < . . . < /i/c(n)( n )- D 

7 Computation of 7r/(x) 

Let / be an arithmetic function, i.e a function defined on positive integers. 
The simplest way to compute tt/(x) defined in (|1.5j) is to generate the primes 
up to x by Eratosthenes's sieve, which is too expansive for large values of x. 

Definition 7.1. An arithmetic function f is said to be completely multi- 
plicative if f(ab) = f(a)f(b) for all a and b. If f ^ 0, this implies f(l) = 1. 

Following ideas of the german astronomer Meissel, Lagarias, Miller and 

/ x 2/s\ 

log a; I 

In this work they also remark that their algorithm allows to compute 7T/(x) 
for every completely multiplicative arithmetic function /. 

This method has been improved in [4J to compute ir(x) with a cost 
/ x 2 / 3 \ 

O 2 — i provided that all the arithmetic operations on integers are of 

\log x J 

constant cost 0(1), not depending on the size of the operands. We show 
here that this improved algorithm may be used to compute vrj(x) whith 

/ x 2/3 \ 

a cost which is still O — ^ — , for a large subset of the set of completely 
\log x J 

multiplicative arithmetic functions. More precisely we have the proposition : 



Odlyzko gave in [8] an algorithm that computes tt(x) with a cost O | 



25 



Proposition 7.1. Let f be a completely multiplicative arithmetic function 
with integer values. Let F be the summatory function of f, 

(7.1) F(a;) = £/(")• 

n<x 

We suppose that all the ordinary arithmetic operations about integers are of 
constant cost 0(1), and that 

1. Each value fin) may be computed in time O(l), not depending of the 
size of n. 

2. There is an algorithm computing 

'X 

, n 



(7.2) S (y,x)= £ /,(n)/(n)F(- 

l<n<y 



in time 0(x 2 / 3 /log 2 x 



Then, there is an algorithm computing it f(x) = > f(p) in time O 



,2/3 



p<.x 
p prime 



log x 



When F(u) can be computed in 0(1) time, the second hypothesis is sat- 
isfied. 

Remarks : 

1. The second hypothesis may seem strange. Let us give a few words of 
explanation. 

• Our computation of vr/(x) begins by choosing y = 0(x l ^ +£ ) . 
Then we compute So = So(y,x) (this is the contribution of ordi- 
nary leaves defined in lemma 5.2, equation (9) in [4] and in lemma 
7.2, equation (7.14) in this article). Function F does not appear 
elsewhere in the algorithm. So being computed, the total cost of 
the other computations is 0[ x 2 ^ 3 / log 2 x ) . Condition (2) ensures 



that our algorithm computes vrj(x) in time 0(x 2 ^ 3 /lo 



,.2 



' X 



• In many cases, F(u) can be computed in time 0(1), then the sum 
(|7.2p can be computed in time 0(y), by precomputing the Mobius 
function, so that the second hypothesis is satisfied. 

2. In Proposition 17.11 we restrict ourseves to the case of integer valued 
functions. The case of real valued functions is more delicate because 
of truncation errors. In [lj, Bach and al. have elaborated an algorithm 
to compute vrj(x) where f(n) = 1/n, and 
x = 1801241484456448 000 = 1.8. . . x 10 18 . 
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Algorithm for 7T/(x) 

We will describe very briefly our algorithm to compute TTf(x), using notations 
and formulas which, when replacing / by 1, reduce to the correponding ones 
contained in [4]. 

For b £ N, let us define b) as the sum of the f(n), for the n's £ [1, x\ 
that subsist after sieving this interval by all primes pi,P2, ■ ■ ■ ,Pbi 

(7.3) 

$(s,0) = f( n ) = F ( x ) and, for 6 >1, $(x,b) = ^ f(n) 

l<n<x i<?i<x 

P-{n)>p b 

For & > 1 and b > 1, let us set 

(7.4) P fc (x,6)= ^ /(n). 

n(n) = fe, F-(n)>p b 

so that, from fjT.3 j) and (I7.4p we get, for x > 1, 

(7.5) *(aj, b) = 1 + Pi (x ; 6) + P 2 (x, &) + ■■■ 
From now on, we choose t/6K 

(7.6) x 1 / 3 < y < \[x and set a = vr(y). 

We will precise later the best choice for y, which is closed to X 1 / 3 . 

Since y > x 1 ^ 3 equation (I7.4p yields P/%(x,a) = for k > 3 and (I7.5P 
becomes 

$(a?,a) = 1 + Pi(x,a) + P 2 (x,a). 
Since Pi(x,a) = ^ /(p) = ^ /(p) = 7T/(x) - 7T/(y), 

(7.7) tt/(x) = *(aj,a)+ir/(y) - 1 -P 2 (x,a). 
Replacing / by 1 (and 7Tj by 7r), formula (|7.7jl is formula (4) in [1]. 

7.1 Initialization of the computation: the 2 basis tables 

After fixing y, by using Eratosthenes's sieve, we precompute the table of 
primes up to y, and the table of the values irf(u) for 1 < u < y. The cost of 
these initializations is 0(ylog 2 y). 
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7.2 Computation of P 2 (x, a) 

Definition (|7.4p and the complete multiplicativity of / give 
P 2 (x,a)= /(M)= E /(?)/(*) 

pq<x PQ^x 

where p and g are primes. The p's figuring in this sum satisfy p < — < — 

Q y 

and we get 

P 2 (x,g)= £ /(p) £ /(<?). 

y<p<x/y p<q<x/p 

X 

We remark that, for p > *Jx, the sum on q vanishes. Since, by (17. 6p . v/x < — , 

y 

we have 

P 2 (x,a)= ^ /(p) ^ /(g) = ^ /(P)(^Q-T/(P-1)) 

y<v<y/x p<q<x/p y<p<Vx 

or 

(7.8) P 2 (x,a)= 53 f(p)n f (y\- 53 /(pK/Cp-l). 

y<p<y f x y<p<V% 

In the above formula, the values of p are bounded above by y/x which is 
larger than y. Thus we cannot find these primes p, nor the values 7Ty(p — 1) 
in the precomputed tables (cf. £17, ip . and we generate them using a sieve 
of [1, \/x\, which we call the auxilliary sieve. The values of x/p lie in the 
interval [1, x/y]. So we will get the values irt(x/p) by an other sieve, the 
main sieve. Let us note that the respective sizes of the sieve intervals, \J~x 
and x/y, are too large to allow a sieve in one pass. Thus the two sieves will 
be done by blocks of size y that must be synchronized. 

• Initialization: Computation of w, the largest prime < \fx and 

of 7Tf(w). By Eratosthenes'sieve we compute the largest prime w 
not exceeding \fx and calculate 7T/(ot). The auxilliary sieve is then 
initialised by putting in the sieve-table the primes p of the block [\fx — 
y + \,i/x[. The main sieve is initialized by sieving the first block 

[A, B] = \y/x, \fx + y — l] . The cost of this phase is O (x 1 / 2 log 2 x^j . 

• Computing P 2 (x,a). We use formula (|7.8p . getting in decreasing or- 
der the primes p €]y,y/x\ and the /(p)tt/(p — 1) from the auxiliary 
sieve, and getting the values 7Ty(x/p) from the main sieve whose suc- 
cessive blocks will cover in ascending order the interval [\/x, x/y). 

We initialize a variable p with the value w, a variable T with the value 
~Kf(w) and a variable P 2 with the value 0. Then, while p > y, we 
repeat : 
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— substract f(p) from T. Thus the new value of T is %f{p — 1). 

— If x/p > B, while x/p > B we replace the block [A, B] by the next 
block [A + y, B + y] and we sieve it. When x/p G LA, -B] we get 
TTf(x/p) in the main sieve table and we add f(p)^f(x/p) — f(p)T 
to P 2 - 

— Using the auxilliary sieve, replace p by its predecessor. 

The final value of the variable P2 is P2(x,a). The first step is negligible in 
cost, compared to the second. Thus the computation of P2(x,a) is of total 

cost 0^1og 2 y^. 

7.3 Computation of $(:r, a) 
The following lemma is proved as lemma 5.1 in [4j. 
Lemma 7.1. For every u > 0, and for b > 1, 
(7.9) $(u,0) = F(u) 



(7.10) *(«, 6) = 6 - 1) - /(p 6 )$ ( -, 6 - 1 

This relation gives an obvious method for computing &(x, a). Starting 
from the tree with the only node <3?(x,a), and applying repeatedly f|T. lOj) we 
get a tree whose all nodes, except the root node, are labelled by a formula 
of the form 

(7.11) M (n)/(n)$g,6) 

where b < a — 1 and n = 1 or n is a squarefree integer with prime factors 

q G {P6+l,---,Po}- 

If we repeat this expansions until all the leaves of the resulting tree are 

labelled by formulas //(n)/(n)$ using (j7.9|) we get the formula : 

(7.12) $(x,a)= M(n)/(n)$(|,0) = £ fi(n)f(n)F g 

l<n<x l<n<x 
P+(n)<y P+(n)<i; 



which, when / = 1 is formula &(x,a) = / i ( n ) ~ ( c f- HI P- 237]) 



1 < n < :c 
P+(n)< H 



The number of terms in (|7.12p is much too large. In order to get a sum 
with fewer terms we replace the trivial rule 

Rule 1 : Expand (f7TTT]) using (f77TU|) if b > 0, 
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which leads to (|7.12p by the new rule 

Rule 2 : Expand node (\7 .lip only if b > and n <y. 

Expanding the computation tree whith rule 2 instead of rule 1 we get 
Lemma 7.2. We have 

(7.13) $(x,b) = S + S, 

where Sq is the contribution of ordinary leaves 



(7.14) 5 = 5 (y,x)= £ M (n)/(n)# (~ Oj = ]T Kn)f(n)F (± 

l<n<y 

anc? 5", the contribution of special leaves, is 



(7.15) S= Kn)f(n)$(-,n(p-(n))-l 



This lemma corresponds to lemma (5.2) in pQ. 

7.3.1 Computation of So 

In the general case, the computation of So is done with a cost 0(x 2 ' 3 / log 2 x) 
thanks to the condition 2 in proposition 17.11 

In the case we will consider later in this work, the computation of i^idix), 

fin) = n, thus F(u) = - - ^~ — - is computed in Oil) time and the compu- 
tation of Sq(x, y) is of cost 0(y) = o(x 2 ^ 3 / log 2 x^j . 

7.3.2 Computation of S 

In the sum (|7.15p , let us set n = rap with p = P~ (n) . Grouping together all 
the n's according to the value of p, we get 

(7.16) S = -J2fiP) E Krn)f(rn)$(^-,n(p)-l)- 

p<y p-(m)> P \mp 



m<y<Cmp 



The computation of S from (|T. 16j) is the complicate part of the algorithm. 
In the following paragraph we show that it is relatively simple to get a cost 
0(x 2 / 3 + e ). 
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7.3.3 How to compute S in (x 2/3+£ 



In this section, we explain a first method to get irf(x), rather simple to 
implement, and whose running time is 0(x 2//3+e ). We take y = x 1 / 3 . Since 
mp > y all the values u = x/mp appearing in (|7,16p are less than x 2 / 3 . We 



sieve the interval 



1, — ) successively by all primes p < y. After the sieve 
V) 

by p, from the definition (I7.3P of <I>, for all the m's such that m < y < mp, 

we get in the sieve table the value <3? ( , ir(p) — 1 ) , and we add to S the 

\mp J 

value /(p)u(m)/(m)5> ( , 7r(p) — 1 

\mp 

But, if we proceed in the naive way, after sieving by each p, we will 
update the sieve table, putting in the case of index u the sum of f(n) for 
the n's, n < u that are still in the table. This is excluded because, for each 
p this would cost 0(x/y) operations, and the total cost of these updatings 
would be 3> ir(y){x/y) = x/logx. As explained in (8] (the 7 last-lines p. 545 
and the first half of p. 546) we use an auxiliary data structure such that, for 
a price of O(logx) time in place of O(l) for each access, we don't need to 
update the sieve table after each sieve. To be a little more precise let us say 
that this structure is a labelled binary tree. There is a leave for each index 
i of the table sieve, this leave is labelled by the value f(i), and each interior 
node is labelled by the sum of labels of its two sons. Proceeding in this way 

the cost of the sieve is O ( — log x log 2 x ) , while the cost of retrieving the 

\V J 

values f(p)fi(m)f(m)& I ,^(p) — 1 ) is O ( it (y)y log x). Both costs are 

\mp J 



0^x 2 / 3+e ^ with our choice y = x 1 / 3 . 
7.3.4 Faster computation of S 



In this section, we explain how to carry out the computation of vrj(x) in 

/ x 2/3 \ 

O I 2 — ] ■ We take y = x 1//3 (logx) 3 log 2 x. To speed up the computation 

\log x J 

of S we partition (|7.16p in 3 subsums S = S\ + S2 + S3, 

Si = - £ M £ mN/M*(^,t(p)-i) 

X a< ~P^y m<y<mp 

.i^A p-(m)> P 



S2 



m<y<.mp 
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S 3 = E MM/(m)*Qjj-,7r(p)-l) 

B <,1 p-( ro )> P vmp 

m<y<mp 

We will show that Si is quickly computed in 0(y) time. S3 will be com- 
puted by sieve, as explained in £17.3.31 but faster because the number of values 
for p is reduced from 7r(y) to ^(x 1 / 4 ). The main part of the computation 
will be the computation of S2. 

As in [4], we first observe that the m's involved in Si and S2 are all prime 
and therefore : 

(7-17) Si = X) /(P) E /(9)*(f .t(p)-i) 

**<P<V P<q - V 

(7.18) S 2 = E /(P) E /(ff)*fe^(P)-l) 

Computing Si As in [4] we remark that, in (|7.17p . we have — < x 1 / 3 < p. 

pq 

Thus, all the values I — ,7r(p) — 1 I are equal to 1. Therefore 

\pq ) 

£1= E /(p)E/(9)= E f(p) My) - */(p)) • 



This value is computed in O(y) additions, using the precomputed table of 
the values Hf{u) for 1 < u < y. 



Computing S3 For each p < x 1 / 4 we precompute the list of all the square- 
free m < y whose least factor is p. 

We sieve the interval 1, ~ successively by all the primes up to x 1 / 4 . 
As soon as we have sieved by p, using the precomputded lists of squarefree 
whose least prime factor is a prime q > p we sum the 



f{p) E i"H/H $ ( — >t(p) - 1 ) 

\ mp J 



P-(m)>p 



for all squarefree m £ [y/p, y) such that P (m) > p. This computation 
is done by blocks, using the auxiliary structure, as explained at the end of 

§E33J 
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Thus the cost of sieving is ( — log x logo x ) . The number of values of p 

\y J 

is -^(x 1 / 4 ) and the number of values of m is less than y, thus the cost of re- 
trieving the values /(p)/i(m)/(m)$ f , 7r(p) — 1^ isO^^x 1 / 4 ) x logx x y^j 

Thus computing S3 is of cost 0\ — log x logo x + yx 1 ^ 4 

\y 

Computing S2 We split the sum (|7.18p in two parts depending on q > 
x/p 2 or q < x/p 2 . It gives 

S 2 = U + V 



with 



and 



u = E /&>) E /(?)*(f.^(p-i) 

1 1 v<q<y \PQ / 



* * P<9<H 

X?<p<x3 q>x /p-2 



v = E /w E /(9)*(f.*(p-i) 

11 P<9 < !/ / 

Computing U With y < ^/x (cf. (|7.6p ). the condition q > x/p 2 implies 
p 2 >x/q>x/y> x 1 / 2 . Thus, 

u= E /(p) E /(9)*feT(p-i) 

From x/p 2 < g we deduce x/pg < p and <fr(x/pq, ir(p) — 1) = 1. and we have 
x/p 2 > p so that 

E E /(<?)= E m(*f<v) 



E7 ^ 



x/y<p<x 1 /' A f-^?7Z,2 \fxJy<p<x x l' A 



Since x/p 2 < q < y the sum {7 is calculated in 0(y) operations with the 
table of values of 717 (w). 

Computing V" For each term involved in V we have p < — < x 1//2 < p 2 . 

Hence, by (|7,3p . <&(x/(pq), 7r(p) — 1) is the sum of f(n) for n satisfying n < 
x/pq and P~(n) > p. These re's are n = 1 and all the primes n satisfying 
p — 1 < n < x/(pq). Thus 

Km' 7 ^ -1 ) = 1 + 7r/ (^) "^(f" 1 ) 
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And we write 
with 



V = V 1 + V 2 



V i = E /(P) E /(9)(1-T/(P-1)) 

x 1 / 4 <p<x 1 / 3 p<q<mm(^,y) 

v* = E m E /(^/(^) 

a; 1 /4<p< a; l/3 p< 9 <min(^,2/) KFH 7 

Computing Vi can be achieved in O(y) time once we have tabulated 7r/(it) 
for u < y. 

Computing V2. We first split V 2 in two parts in order to simplify the 
condition q < min(x / p 2 , y) : 



v*= E /(p) E (^)+ E /(p) E /(*)*/( 

^ 1/4 <P<^f P<9 -^ J%<P<^ 3 P«l<fl 



X 

pq 



In the purpose to speed up the computation of the above two sums we now 
write, V 2 = Wi + W 2 + W 3 + W A + W 5 with 



wi = E /(p) E (i) 

sr 

E /(p) E /(*)*/ 

E /(p) E 

^ = E /(p) E ^^(^ 
E /(p) E 



1/4 <p<^ p<q<v 

Wo = 



W- 



<p<x!/3 ,/-<<?< — 



Computing Wi and W2 These two quantities need values of Wf(x/pq) 
with x 1 / 3 < x/pq < x 1 / 2 . These are computed with a sieve of the inter- 
val [1, y/x\. The sieving is done by blocks, and, for each block, we sum 
f{p)f(q)^f(x/pq) for the pairs (p, q) satisfying the conditions of the sum W\ 
or W 2 and such that x/pq lies in the block. 
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The cost of this computation is the sum of three terms 



The cost of the above sieve on [l, y/xj is Oi^/x log 2 x) . 



• The cost of adding the terms of the sum Wi, O 

• The cost of adding the terms of the sum W2, O 



y log 2 x 

x 3 / 4 \ 
yV4 log 2 x I ' 



Computing W3 and W5 For W3, for each p we apply lemma [2761 with z 



x/p and u = y. Thus, for each value of p, the sum on q costs Oyir^x/p) 
and the total cost of the computation of W3 is 







( 

E - 



For W5, for each p we apply lemma ETUI with z = x/p and u = x/p 2 . Thus, 



for each value of p, the sum on q costs Oyir(y^x/p)j , and the total cost of 
the computation of W5 is 







E 



7T 



i<p<I 1 /3 



Thus the costs of computing W3 and W§ add to 








E 



log - 

y^ <p < x i/3 \ /' 



x 2/3 

log 2 X 



Computing VF4 We simply sum over (p, q). There would be no advantage 
to proceed as for W3 since most of the values ftf{x/pq) are distinct. The 
cost is 

/ \ 



O 



E 



7T 



01 



2/3 



log x 



As in [3], section 8, we then see that, since y = x 1 / 3 log 3 x log 2 x, the 
total cost of the computation of vr/(x) is 0(x 2 / 3 / log 2 x) . 
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8 The algorithm to calculate h(n) 
8.1 The function G(pk,m) 

The function G(pk,m) has been introduced and studied in [5]. 

Definition 8.1. Let pk be the k-th prime, for some k > 3 and m an integer 
satisfying < m < Pk+i — 3. We define 

(a l\ ni \ Q1Q2 ■ ■ ■ Qs 

(o.l ) G{pk,m) = max ■ 



qiQ2 ---qs 



where the maximum is taken over the primes Q±, Q2, . . . , Q s , qi, q2, . ■ ■ , q s 
(s >0) satisfying 

(8.2) 3 < q s < q s -i < . . . < q x < p k < p k+ i < Qi < Q2 < ■ ■ ■ < Qs 
and 

s 

(8.3) J^(Q 4 - Qi ) < m. 

i=l 

The additive function t (cf. §1. ip can easily be extended to fractions by 
setting 

£(M/N) = £(M) - £{N) 

when M and are coprime or are both squarefree. Therefore, the inequality 
(|8.3p implies 

(8.4) i(G( Pk ,m))<m. 
8.1.1 Properties of G(p k , m) 

Obviously, G{p k ,m) is non-decreasing on m and G(p k , 2m + 1) = G(p k , 2m). 
The maximum in (|8,ip is unique (from the unicity of the standard factoriza- 
tion into primes). For small m's, we have 

(8.5) 0<m<p k+1 -p k G(p k ,m) = l. 
From Proposition 8 of [5j, we have 

/o c ^ Pk+1 , n ; n - Pk+1 

8.6) < G(p k , m) < . 

(Pk+i - m)* p k +i - m 

Note that if p k +i — m is prime, then (j8.6|) yields the exact value of G(pk,m). 
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8.1.2 Computation of G(pk, m) 



In [5j §8], two algorithms are given to calculate G(pk,m). 

The first one is a combinatorial algorithm. In its first step, the primes 
allowed to divide the denominator of G(pk,m) are determined. From (I8.2|) 
and (|8.3p . they are all the primes in the range [(p^+i — m), pk], say P\ < 
Pi < ... < Pk- Similarly, the primes authorized to divide the numerator 
are all the primes Pr+i < Pr+2 < ■ ■ ■ < Pr in \pk+i, Pk + m]- By setting 
p' = {P 1 ,P 2 , . . . , P R }, from the definition (fL14")) . we get 

G{Pk,m) = h K {P 1 +P 2 + ... + P K + m,P') 

"l "2 ■ ■ ■ tk 

and hj(n,V) can be computed by induction on j in a way similar to that 
exposed in §1.41 In [5J §8], one can find the details and also some tricks to 
improve the running time of this combinatorial algorithm which, however, 
remains rather slow when m is large. 

The second algorithm, which is more sophisticated, is based on the fol- 
lowing remark : if G(pk,m) = Q 1 ^ 2 ^£ anc l m [ s large, the least prime 

QiQ2 • • • q s 

factor q s of the denominator is close to Pk+i — m while all the other primes 
Qi,..., Q s , qi, ■ ■ ■ , q s -i are close to p k . 

More precisely, the following proposition (which is Proposition 10 of |5j) 
says that if pt+i — m + 5 is prime for some small 5 and if G(pk+i, S) is not 
too small, then the computation of G(pk,m) is reduced to the computation 
of G(pk+i,m') for few small values of mf , which can be done by the above 
combinatorial algorithm. 

Proposition 8.1. We want to compute G(pk,m) as defined in (|8.ip with p^ 
odd and pk+i — Pk < m < Pk+i — 3. We assume that we know some even 
non-negative integer 5 satisfying 

(8.7) Pk+i — m + 5 is prime, 



G( Pk+1 ,5)>l+ d 



Pk+i 
and 

(8.9) s < £ < Ban. 

If 5 = 0, we know from (18.61) that G(pk, m) = — Ehll If 6 > 0, we have 

Pk+i - m 

(8.10) G(p k ,m) = max G(p k+ i, m - Pk+i + q), 

q prime q 
Pk+l-m <q<q 
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where q is defined by 

/ ei -n - Pk+iPk+2(Pk+i -m + S) 35 

(pjfc+i + d)(pfc+i - 35/2) 2 

How to compute G(pk,m) ? The combinatorial algorithm should be tried 
if m is small, but it is quadratic in m and has no chance to terminate if m 
is larger than, say, 10 6 . We have no guarantee that the conditions of Prop. 
I8.ll are satisfied. However in all our numerical applications, we have found 
5 < 1000 in ||g2J (see [SJ §9.2]), so that, by fl8~T0]) and (fSTTTj) . we have 

"i- - Pfc+i + g < m - pfc+i + 9 < Pfc+2 - Pfc+i + y 

and, in (|8.10p . G[pk+l, m — Pk+i + (?) can be easily calculated by the combi- 
natorial algorithm. 

8.2 Description of the algorithm to compute h(n) 

To compute h(n), the first step is to determine pk and Ok defined by (|1.8p . 
This step is explained in § 18.2.11 and will furnish also Pk+i and n' = n — pk- 

8.2.1 Computation of pk and Ok 



1. Compute x = y Li 1 (n), so that Li(x 2 ) = n and x ~ \fn logn. 

2. Using Prop. 17.11 we compute iiid{x) in time 
o(x 2 / 3 /log 2 x) = oU/z/ilogn)- 5 / 3 ). 

3. To get <7fc, we have to add (if iTid(x) < n) or to subtract (if TTid(x) > 
n) to Hid{x) the primes between x and calculated by sieving. In 
practice, this step is very short. But we are able to estimate it only 
under Riemann's hypothesis. By lemma [231 we have 



Li(pl)-^pl /2 log p k <cr k <n< cT fc+ i < Li(p| +1 )+— pj£ logpfc+i 
which implies 



2 

n = U(jpl) + O (pf 2 logpk) ~ Li(p 2 ^ - '" r 



21ogp fc 

Therefore, we get logn ~ 2\og pk, Pk ~ y/nTogn and 
(8.12) Li(^) =n + o(n^ 4 (logn) 7 ^ 



Further, since x ~ pk ~ \/n logn, we have 



ln-Li(p 2 )| = |Li(x 2 )-Li(p 2 )U ^— ^ ~2|x-p fc ' ; 



so 



21ogx Y logn' 

that, from (|8~T2l) . |vr(x) - 7r(p fc )| < \x - p k \ = C>(n 1/4 (logn) 9/4 ) . 
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8.2.2 Computation of h(n)/N k 



By Corollary 16.11 we have h(n) = h k (n). Let us set n' = n — a k . If n' = 
Pk+i — 1 or n' = Pk+i — 2, Proposition 15. 31 yields h(n) = N k+ \/2. So, we may 
suppose n' < <r k+ i — 3. From the definition (|8.1|) of function G, we have 

(8.13) h(n) = h k {n) = N k G( Pk , n') 

and we compute G{p k ,n') as explained in £18.1.21 In practice, the compu- 
tation of G(p k ,n) is very fast. However, as explained in §8.1.21 we have no 
estimation of the running time. 

h(n) 

Below, are listed some values of = G(p k ,n') together with p k , 

N k 

n' = n — a k , e = e(n) the largest integer such that h(n — e) = h{n)) and, 
if the algorithm of Proposition 18.11 is used, 5 and Q, the number of primes 
used in the sum (|8.10p . 

n = 10 12 , p k = 5477081, n = 4935150, e = 0, 5 = 18, Q = 1, 
29998525822277 5477089 x 5477093 



G(p k ,n') 



2968309525031 5477081 x 541951 



n = 10 35 , p k = 2898434150644708999, n = 1886081812111845520, e = 16 

. , nA „ r , N 2898434150644709023 

5 = 134, Q = 5, G(pu, n ) = 

' ' ' 1012352338532863519 

The values of /i(10 a ) for a < 35 and of h = 2 b for b < 116 can be found 
on the authors's web sites [21 ED]) together with the Maple or Sage programs 
computing h(n)/N k . 



9 An open question 

Given n and j < k(n), how to compute hj{n)l We have not succeeded in 
solving this problem when n is too large to use the naive algorithm described 
in 11.41 The case j = 2 is already not that simple. 

A first step is certainly to calculate r = r(n,j) defined by (|4.ip . which 
can be done by the method of § 18.2.11 If we are lucky enough that q = 
Pj+ r +i—n' = pj +r+ i — (n — (jj +r + a r ) is prime, then the value hj{n) = jV ^' +1 
is given by (|4.24j) . 

In the general case, by setting n' = n — (cJj+ r — o>), one may think 
that h'j(n) = G(pj+ r , n) has a good chance to be the value of hj(n). 

However, there are exceptions. 
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J = 


1 


2 


3 


4 


5 


6 


n = 2 


2 












3 


3 












4 


3 












5 


5 


6 










6 


5 


6 










7 

8 


7 
7 


10 
15 










9 


7 


15 










10 


7 


21 


30 








11 


11 


21 


30 








12 


11 


35 


42 








13 


13 


35 


42 








14 


13 


35 


70 








15 


13 


35 


105 








16 


13 


55 


105 








17 


17 


55 


105 


210 






18 


17 


77 


110 


210 






19 


19 


77 


165 


210 






20 


19 


91 


165 


210 






21 


19 


91 


231 


330 






22 


19 


91 


231 


330 






23 


23 


91 


385 


462 






24 


23 


143 


385 


462 






25 


23 


143 


455 


770 






26 


23 


143 


455 


1155 






27 


23 


143 


455 


1155 






28 


23 


187 


455 


1365 


2310 




29 


29 


187 


715 


1365 


2310 




30 


29 


221 


715 


1365 


2730 




31 


31 


221 


1001 


1430 


2730 




32 


31 


247 


1001 


2145 


2730 




33 


31 


247 


1001 


2145 


2730 




34 


31 


253 


1001 


3003 


4290 




35 


31 


253 


1309 


3003 


4290 




36 


31 


323 


1309 


5005 


6006 




37 


37 


323 


1547 


5005 


6006 




38 


37 


323 


1547 


5005 


10010 




39 


37 


323 


1729 


5005 


15015 




40 


37 


391 


1729 


6545 


15015 




41 


41 


391 


2431 


6545 


15015 


30030 


42 


41 


437 


2431 


7735 


15015 


30030 


43 


43 


437 


2717 


7735 


19635 


30030 


44 


43 


437 


2717 


8645 


19635 


30030 


45 


43 


437 


2717 


8645 


23205 


39270 


46 


43 


493 


2717 


12155 


23205 


39270 


47 


47 


493 


3553 


12155 


25935 


46410 


48 


47 


551 


3553 


17017 


25935 


46410 


49 


47 


551 


4199 


17017 


36465 


51870 


50 


47 


589 


4199 


19019 


36465 


51870 



Figure 1: Table of hj(n). 
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